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Abstract 

Resonance searches generally focus on narrow states that would produce a sharp peak 
rising over background. Early LHC running will, however, be sensitive primarily to broad 
resonances. In this paper we demonstrate that statistical methods should suffice to find 
broad resonances and distinguish them from both background and contact interactions 
over a large range of previously unexplored parameter space. We furthermore introduce 
an angular measure we call ellipticity, which measures how forward (or backward) the 
muon is in rj, and allows for discrimination between models with different parity vio- 
lation early in the LHC running. We contrast this with existing angular observables 
and demonstrate that ellipticity is superior for discrimination based on parity violation, 
while others are better at spin determination. 



1 Introduction 



With the LHC running smoothly and the "re- discovery" of the Standard Model (SM) well 
underway, it is time to contemplate what new physics this early LHC run might access. A 
significant amount of data (~ 1 fb~ ) will be accumulated at y/s = 7 TeV, and it is important 
to consider what kinds of physics beyond the SM can be studied at this energy and luminosity, 
particularly with several years' wait until the anticipated 14 TeV run. The models of greatest 
interest in this regard are those that will be accessible to early LHC running but are outside 
of the range of the Tevatron. 

The first and simplest target for new physics searches will be resonance searches. The 
Drell-Yan channel into leptons will likely be the simplest channel to study, since the SM 
contributions are well-understood, the background is very low at high invariant mass, and 
many models of new physics have sizeable branching fractions to leptons, particularly muons. 
Our focus will therefore be on muons, though we will also briefiy consider the additional 
contributions from electrons or photons that could be present and increase discovery reach. 
Muon detector efficiencies and resolutions are very good at both CMS and ATLAS, as are 
those for electrons. However, models involving strong electron couplings are generally tightly 
constrained by LEP, making dimuon searches more interesting in those cases. The production 
of strongly-interacting final states (such as tt) could also be of interest, as they are enhanced 
by a color factor, but this advantage is typically negated by higher QCD backgrounds, lower 
efficiencies, and worse resolution, so we will not consider these further. 

The pp nature of the Tevatron ensures that it will serve as the best probe of neutral res- 
onances below 1 TeV. The LHC will nonetheless dominate for other cases due to the higher 
center-of-mass (CM) energy. This makes high-mass resonances, which were out of the Teva- 
tron's reach, kinematically accessible. Furthermore, the large gluonic parton fraction at low 
X means that resonances coupling to gluons (such as Kaluza-Klein (KK) gravitons in the 
Randall- Sundrum (RS) framework) will have enhanced cross sections over the Tevatron. Even 
so, to be visible with only 1 fb~^ of data, large couplings to leptonic final states are most likely 
essential. Since the width of a resonance scales as the square of the coupling, this implies that, 
during the early LHC running, we will find broad resonances with large decay widths. 

Broad resonances are more challenging to study than narrow ones. With narrow reso- 
nances, the invariant mass of the signal events are tightly clustered around the particle's 
mass. For muons, the detector resolution is sufficiently good that performing a simple count- 
ing experiment in the region around the resonance pole (with a width on the order of the 
muon resolution) often gives a signal-to-background ratio that is high enough for detection. 
Detecting broad resonances will be more challenging since signal events are spread out over 
a wider area so the signal will not be as clearly identifiable by its shape. Distinguishing sig- 
nal from background is further complicated by parton distribution uncertainties in the high 
invariant mass region, as well as uncertainties in the integrated luminosity, both of which are 
often estimated using the Drell-Yan process that we're using as a probe of new physics. 

The question then becomes whether or not we can distinguish these events from background 
or other forms of strongly-interacting physics such as contact interactions. In this paper we 
show that characteristic features of the invariant mass distribution will be sufficient to find 
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and identify resonances, even at the first run of the LHC In its simplest form, our method 
involves looking for an absolute rise in the differential cross section, or "upturn" , which would 
clearly distinguish any new physics from the SM or a contact interaction, both of which predict 
a falling distribution. Our fuller, more sophisticated statistical analysis uses the maximum 
likelihood method to distinguish between resonances and contact interactions (of which the SM 
is the limiting case where we take the contact interaction scale A -> oo) and covers a greater 
parameter range. We demonstrate that differentiating a broad resonance from background 
and contact interactions is possible over much of the accessible range of parameter space with 
a significant number of signal events. 

We also use angular distributions to further distinguish new physics models from one an- 
other and from the SM. We study the muon pseudorapidity (r/") distribution, and demonstrate 
that a variable we call ellipticity can be used to distinguish among different chiral structures. 
We find that ellipticity is good at distinguishing models with different parity violation, whereas 
observables based on the muon angle in the CM frame (6*) are more useful for spin deter- 
mination. Model discrimination based on angular information generally requires more events 
than will be accessible in the first LHC run for resonance masses larger than 1 TeV, so we also 
consider the prospects of model discrimination with y/s = 10 TeV and different integrated 
luminosities. 

Our analysis follows others who have investigated the LHC reach of Z' gauge bosons during 
the a/s = 7 TeV run with various integrated luminosities [H HI [3l HI [5]. We extend these 
analyses to include searches for RS KK gravitons, and widen the study to include strongly- 
coupled, broad resonances. 

2 Resonance or Contact Interaction? 

We begin our work by motivating the fundamental question addressed in this paper: can we 
distinguish resonances from contact interactions in the 7 TeV run of the LHC? For concrete- 
ness, we examine Z' models and KK gravitons within the RS framework (we will describe the 
details of these models in section [3]). As we will demonstrate in section HI the largest accessible 
regions of parameter space in the early running of the LHC will be strongly coupled, resulting 
in a broad resonance. 

As the coupling of the theory increases, the resulting resonance broadens and flattens, 
losing its distinctive "bump", and instead begins to look more and more like background, a 
contact interaction, or some other new physics that predicts excess events at high invariant 
mass. Nonetheless, discriminating a resonance from other possible models (including the SM) 
should be possible for many of the broad resonance models we study, although it is not a 
priori obvious that this is the case. To show why model discrimination may be difficult, we 
plot the differential cross section for both RS and Z' resonances in Fig. [H varying the mass 
and coupling to get a qualitative sense of how the invariant mass distributions change with 
these parameters. We also plot in Fig. |2] the differential cross section for a set of contact 
interactions for comparison with the resonances. 

To illustrate the aforementioned difficulties in distinguishing different models of new physics 
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Figure 1: Differential cross sections for RS (top) and non-universal Z' (bottom), as defined 
in section [3l Top left: graviton of mass 2 TeV with k/Mpi = 0.7 (highest), 0.5, 0.3, and 0.1 
(lowest). Top right: k/Mpi = 0.4 and masses 1.4 TeV (left), 1.9 TeV, 2.4 TeV, and 2.9 TeV 
(right). Bottom left: Z' of mass 2 TeV with e = 1.4 (highest), 1.1, 0.8, and 0.5 (lowest). 
Bottom right: Z' with e = 1 and masses 1.4 TeV (left), 1.9 TeV, 2.4 TeV, and 2.9 TeV (right). 
The SM is shown in black on all plots. 
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Figure 2: Differential cross sections for four different contact interactions at ^/s = 7 TeV. The 
models are given by Eq. f|T71) with various signs of t]. From highest to lowest: VV constructive 
(green), LL constructive (blue), LL destructive (orange), VV destructive (purple), SM (black). 
The contact scale is A = 4 TeV for the VV models and A = 2.87 TeV for the LL models. The 
energy scales are chosen to have the same high-energy behavior for all four. 
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Figure 3: Illustration of the potential difficulty of determining new model type from data, 
whether it be a resonance or a contact interaction. Data was simulated for RS (left) and non- 
universal Z' (right), and is shown in blocked red. Best fits to data from resonance and contact 
models are shown in red and blue dots, respectively. Black dots show the SM prediction. 



through the invariant mass distribution, we ran two simulations showing possible outcomes 
of an experiment at y/s = 7 TeV and 1 fb~^ (see Fig. |3]). We fit a resonance and a contact 
interaction to the data and, by sight, it seems that both might be consistent with the data up 
to statistical fluctuations. We want to see when it is possible to reliably distinguish between 
the two types of new physics, and how this discriminating capability varies with the parameters 
of the theoryll] In section El we will look at two methods for doing so. The ffist is relatively 
simple, looking for an absolute rise in the differential cross section. The second employs a 
statistic that is commonly used in experiments with few data points in each bin. We show 
that both of these methods can be used to distinguish resonances from contact interactions 
for most models predicting > 5 events in the early LHC running, though the more careful 
statistical method covers a larger range of parameters. 



3 Models of new physics 

In this section we define the two types of resonance models we use to demonstrate our methods, 
namely Z' models of a new U{iy gauge boson, as well as RS models with at least one accessible 
KK graviton resonance. We also define the contact interaction models that we will compare 
with the resonances. 

^If we take the resonance mass and coupling to be very large, we can integrate out the new physics and 
we do end up with a contact interaction. This, however, is not what we're talking about here; rather, we are 
concerned about the case where the observed signal events are in the region around the resonance peak. 
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3.1 Z' models 

Z' bosons commonly appear in physics beyond the SM (for a recent review, see [6]). The 
couphng of a Z' to the SM can be written as 

liz>=gz'Z'^rz: (1) 

where 

J'z' = Y.^{f)Hf. (2) 

/ 

The charges will be discussed more in-depth shortly, although following j^, we will focus 
on Z' that couple to B — L and B — 3L^. We define the gz' coupling in terms of the gz 
coupling by introducing the parameter e: 

gz' 

e = — . 3 

gz 

We compute the differential and total cross sections to determine how these quantities 
scale with the model parameters. The partonic cross section has the form 

^ ~ ^^T-^ . ^9 xn ; — ^ „ + interference, (4) 

(s - M|,)2 + Imn(s)2 ^ ^ 

where near the resonance or in the narrow width approximation, we have a partial width to 
states of mass m 

r ~ Mz' 



l + O 

The partonic cross section evaluated at the pole is 



2 

m 



AH, 



(5) 



ml) - 4, (6) 

and so the peak cross section decreases with higher M^/, as expected. In a hadron collider, 
however, the partonic cross section, a{s]ij — )■ = a'ij{s), is integrated against the corre- 

sponding parton luminosity function, dCij/dr, defined as follows: 

dr, 1 dr 

-^(r) = — ^ / -[/.(x)/,(r/a;) + /.(r/x)/,(x)], (7) 
rfr 1 + Oij X 

where r = s/s, to give the hadronic cross section 



^hadronic ( ^\ \ ^ / U 

dr 



' . dCi 



The hadronic cross section can be evaluated analytically only in the narrow width approxi- 
mation, giving for CM energy ^/s 

^ ^ — ■ -ri^z'), (9) 

s dT 




Figure 4: Parton luminosity dCu/d'T (t 
the LHC, = 7 TeV. 



s/s) for uu contribution to scattering processes at 



where tz' = M^,/s. This is independent of the mass Mz', except for the dependence on tz' in 
the parton luminosity function. The parton luminosity function decreases monotonically with 
r, as seen in Fig. Ho 

Z' phenomenology is dictated by the parameters Mz' and e. In our analysis, we consider 
all values of these parameters consistent with experimental constraints. To observe higher- 
mass resonances, the coupling e must be sufficiently large to compensate for the small parton 
luminosity, meaning that the Z' models most visible during the early LHC running will be 
broadjfl 

The only rigorous constraint on Z' couplings comes from anomaly cancelation requirements. 
These can be satisfied in myriad ways - for instance, by the presence of hidden exotic states 
- but the simplest are "minimal models", which consist purely of couplings to the SM plus 
right-handed neutrinos [3], [5]. The minimal models allow us to adequately study the various 
Z' possibilities in a systematic way without introducing too many additional parameters, and 
we will restrict our attention to these. Anomaly cancelation then limits the overall charges 
to a linear combination of two independent charge assignments: hypercharge and B — L. 
The hypercharge component leads to mixing with the Z and is tightly constrained by LEP 
measurements. In order to study Z' models with stronger coupling to the SM, we will focus on 
Z' of the B — L type. Strictly speaking, the couplings will mix under renormalization group 
running [3], so models that unify with the SM at the GUT scale will generically see both 
types of charges, but this effect is small and, for the purposes of LHC phenomenology, we can 
ignore higher energies and simply examine the couplings at the TeV scale. 

^We probe the regime where the luminosity can no longer be taken as a simple power law: the total partonic 
luminosity from uu collisions at the LHC is approximately given by the expression dC/dr — 0.036t^-^ °^(1 — 
T)0-02exp(-27.2T). 

^Strictly speaking, this scaling argument only applies in the narrow width approximation, but wc find 
a similar approximate power law scaling of hadronic cross section with e even at large coupling, and our 
conclusion remains the same. 
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Although not as tightly constrained as a Z' that couples to hypercharge, B — L models are 
also strongly constrained by electroweak data if they are flavor-universal due to the constraints 
on electron-quark 4-fermion operators from LEP. If we only consider flavor-universal models, 
the LEP constraints are sufficiently strong to exclude all broad resonances with (9(TeV) mass, 
and since these are the only type that will be visible in early LHC running, there will be no 
new parameter space open for exploration. A way around this is to couple the Z' exclusively 
to muons in the leptonic sector [H [7j; the theory will still be anomaly-free if the charges 
are B — 3L^. These non- universal models evade the bounds from LEP while giving large 
branching fraction to muon final states, enhancing the visible signal at the LHC. This means 
that there are regions of the parameter space that are accessible at the early LHC but that 
have not already been excluded at LEP. Our analysis generalizes that of [1] pertaining to Z' 
coupling to B—3L^, considering in particular the case of broad resonances and its experimental 
implications. 

3.2 Brane RS models 

Many of the best-motivated theories of new physics are those providing a solution to the 
hierarchy problem. One class of examples are RS models, in which the hierarchy between the 
Planck and electroweak scales is established with the addition of a warped extra dimension 
and the Higgs field confined to a brane at the IR scale [5]. For theories with the entire SM 
confined to the TeV brane, KK modes of the graviton could be detectable in early running if 
the TeV-brane energy scale is sufficiently low. 

The couplings of higher KK modes of the graviton are suppressed by the warped scale 
A,r = Mp\e~^''^''. Here, k is the AdS curvature scale and is the radius of the compact 
dimension. If the hierarchy problem is solved, then A^ ~ TeV. The Lagrangian for the 
coupling of the KK graviton modes is 



with energy-momentum tensor T^"^ and n-th KK graviton hl^J . The free parameters in this 
model can be taken to be k and A^^; for phenomenological reasons, it is more convenient to 
exchange these for the dimensionless parameter k/Mp\ and the first KK graviton mass 



where Xi = 3.83171 is the first zero of the Bessel function Ji{x). The subscript g is used 
for the first mode because it is the graviton mode that will first appear at the LHC, and 
higher modes will likely be too heavy for discovery in the early running. These parameters 
are better-motivated from the standpoint of a resonance search and are directly analogous to 
the parameters e and Mz' in Z' models. 

As with the Z', we compute the partonic and hadronic cross sections for lepton production- 
here via KK graviton exchange. For both gluon- and quark-initiated processes, the partonic 




(10) 




(11) 
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cross section takes the form 



1 

■ (s-M2)2 + Imn(s)2' ^^^^ 

where n(s) is the graviton self-energy. On-resonance, or in the narrow width approximation, 
the self-energy contribution can be simplified to give 

Im n(s) ^ Im n(M|) = Mg T, (13) 

where the graviton partial width to final states of mass m scales as 



k 



The partonic cross section evaluated at the pole is 



(14) 



^(K) - -^^^ (15) 

and the hadronic cross section is 

a -(rg), (16) 

where Tg = M^/s. This behavior is clearly analogous to that for a Z', and we are likewise led 
to the conclusion that strongly-coupled physics and broad resonances will appear first at the 
LHC. _ 

Many analyses, beginning with [11], consider only RS models with coupling k/Mpi < 
0.1 due to perturbativity constraints, restricting their studies to narrow resonances. This 
constraint, which comes from bounding the tree-level AdS curvature |i?5| = 20k^ < M^, 
neglects loop effects and is likely overly conservative. 



3.3 Contact interactions 



We now discuss contact interactions in more detail, as we saw in section [2] that, with small 
statistics, they might be confused with broad resonances. Contact interactions are convenient 
for parameterizing the low-energy behavior of new physics contributions. We consider generic 
contact interactions that are consistent with current experimental data and remain agnostic 
about possible UV completions of the theory. 

The contact interactions most relevant to us are those involving quarks and muons. They 
are defined by a scale A, and we use the parameterization [12] 



An 
A2 



+ '7RL (qRi'^qR) (P'LIuIJ'l) + vrr (qR^qR) {^rIuIj^r) 



(17) 
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where |?7| = 1 or and can give chiral or axial/vector-like couplings. 

Models with a vector-like coupling, denoted by VV, will have all couplings set to 1 or — 1, 
each with the same sign. A left-handed model, denoted by LL, has only one non-zero coupling, 
?7ll, and thus only couples to left (right) handed quarks (anti-quarks) and muons (anti-muons). 
For these two models, the signs of the ?7's determine the nature of the interference with the 
SM. The SM is recovered in the limit A — > oo. 

The contact interaction ( 1T71) arises from integrating out some spin-1 degrees of freedom. 
This seems like the best choice for mimicking a Z' resonance, but at first glance, it may seem 
that a contact interaction generated by integrating out spin-2 degrees of freedom might do 
better at faking a graviton resonance. Such contact interactions, however, are dimension 8 and 
suppressed by extra powers of the scale A. As a result, the event rate from such operators will 
not compete with the event rate from a broad resonance without running into serious conflict 
with existing experimental bounds. Thus, we find that even for RS gravitons, the interaction 
f lT7|) is the best basis for comparison (we may also wish to consider scalar-type couplings, but 
as in the high energy limit, this only affects the angular distribution, we do not consider them 
here) . 

4 LHC dimuon cross sections 

The current goal for early LHC running is 1 fb~^ at = 7 TeV. To get a sense of the ultimate 
LHC reach, however, and also to see qualitatively what the possible outcomes of resonance 
searches may be at later times, we consider both the early run and also the cross sections at 
\fs = 10 TeV and ^/s = 14 TeV, even though the analyses presented in this paper are mostly 
concerned with the first case. We consider the RS and Z' models discussed in section [3] for 
different values of the mass and coupling. As anticipated in section [21 we demonstrate that 
the largest cross sections are in regions of parameter space with broad resonances, leading to 
the question of model discrimination between the broad resonance and a contact interaction 
that will be addressed in section [51 

We present our results for Z' signal cross sections in Fig. [51 and [6] (for the flavor universal 
and non-universal cases, respectively), and our results for RS signal cross sections in Fig. [71 
We mark contours of constant cross section, including the contours that will give 5 events 
at certain benchmark luminosities: 1 fb~^ at 7 TeV, 10 fb~^ at 10 TeV, and 100 fb~^ at 14 
TeV. We indicate constraints using dashed lines. Details of the calculation and the relevant 
direct and indirect constraints follow later in this section. In the plots shown here, the regions 
of parameter space most ripe for exploration in the Drell-Yan channel at the LHC are those 
that couple strongly to muons, giving rise to broad resonances. The only exceptions are the 
Z' models that are pure B — L: since they couple universally to muons and electrons, the 
broad resonance regions are largely excluded by LEP. Indeed, even at high CM energy and 
luminosity, a. B — L Z' will appear as a narrow resonance and the standard search strategies 
will apply. 

We compute the differential cross sections at leading order in a using the MSTW 2008 
PDFs [13], with invariant-mass-dependent K- factors applied to quantify NLO QCD effects 
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Figure 5: Plots showing cross sections to dimuon final states for a Z' with B — L charges to 
SM. The cross sections are shown for -y/i = 7 TeV (above left), 10 TeV (above center), and 14 
TeV (below). The solid line cross section contours are, from top to bottom: 1 pb (green), 100 
fb, 10 fb, 1 fb, and 100 ab (brown). Dashed curves indicate the cross section for 5 events at 
certain benchmark luminosities: 1 fb~^ at 7 TeV, 10 fb~^ at 10 TeV, and 100 fb~^ at 14 TeV. 
The shaded region is excluded by LEP and the Tevatron. 



[m [15]. Other necessary cuts due to detector geometry or triggering, for example \ri\ < 2.5 
and pt > 20 GeV, are incorporated into the computation by employing an invariant-mass- 
dependent acceptance, following [3l[T6]. The total cross section is then found by integrating 
the invariant mass distribution in some window. Typically, the width of this window is some 
multiple of the resonance width F, but when we look at very broad resonances with M ~ F, 
this will include events at very low invariant mass where SM background dominates. We 
therefore instead choose a lower cut-off for our window of Max(M — 2F,M^,), where is 
the invariant mass at which the signal differential cross section is twice as large than the 
background. Our upper cut-off for the window of integration will be ^/s, the CM energy of 
the collider. 

We do not exclude any regions of parameter space solely for reasons of theoretical prejudice, 
but consider only experimental constraints. The most important ones for our study will be 
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Figure 6: Plots showing cross sections to dimuon final states for a Z' with B — 3L^ charges 
to SM. The cross sections are shown for ^/s = 7 TeV (above left), 10 TeV (above center), and 
14 TeV (below). The solid line cross section contours are, from top to bottom: 1 pb (green), 
100 fb, 10 fb, 1 fb, and 100 ab (brown). Dashed curves indicate the cross section for 5 events 
at certain benchmark luminosities: 1 fb~^ at 7 TeV, 10 fb~^ at 10 TeV, and 100 fb~^ at 14 
TeV. The shaded region is excluded by LEP and the Tevatron. 



the bounds on 4-fermion operators from LEP [T7] and the Tevatron [18]. The bounds from 
LEP in particular strongly constrain new physics coupling to electrons, such as the Z' coupling 
to B — L. One exception is the RS KK graviton that we have been discussing. It couples 
to the energy-momentum tensor and the lowest- dimension operator contributing to s-channel 
graviton exchange is dimension 8, which is much less strongly constrained than the dimension-6 
operators generated by massive vector exchange. 

We also consider the strongest direct bounds on neutral resonances, which come from 
the Tevatron [191 120]. These direct bounds do not extend much beyond 1 TeV because the 
Tevatron parton luminosity falls off very rapidly above this scale. Finally, we incoporate the 
bounds from muon g — 2 experiments [211 ESI EH] , although these do not strongly constrain 
RS models due to the cut-off dependence of the graviton contribution to g — 2, and so there 
are many values of Mg and k/Mpi that can give the correct contribution. For Z' models, the 
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Figure 7: Plots showing cross sections to dimuon final states forfor the first KK mode of the 
graviton in RS models. The cross sections are shown for ^/s = 7 TeV (above left), 10 TeV 
(above center), and 14 TeV (below). The solid line cross section contours are, from top to 
bottom: 1 pb (green), 100 fb, 10 fb, 1 fb, and 100 ab (brown). Dashed curves indicate the 
cross section for 5 events at certain benchmark luminosities: 1 fb~^ at 7 TeV, 10 fb~^ at 10 
TeV, and 100 fb"^ at 14 TeV. The shaded region is excluded by LEP and the Tevatron. 



constraints from 4-fermion operators are stronger than those from g — 2. 



5 Shape discrimination 
5.1 Upturn analysis 

As a simple first attempt at studying the shape of the distribution, we take advantage of one 
of the defining characteristics of a resonance: namely, the accompanying "bump", or absolute 
rise and fall in the cross section. The background (or a contact interaction) will never actually 
give an increase in the cross section. We therefore look at regions of parameter space where 
we can, with 5a confidence, say that there is a rise in the cross section over a flat background, 
indicative of a resonance with reasonable certainty. We find that for both KK gravitons and 
Z' , the upturn analysis works for masses up to 1.5-2 TeV, and widths up to about 40% for Z' 
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Figure 8: Regions of parameter space for non-universal Z' (left) and RS (right) models acces- 
sible (5cr) at the early LHC via the upturn analysis are shaded in blue. Cross-section contours 
are shown in dashed black lines: 1 pb (top), 100 fb, 10 fb, 5 fb, 1 fb, and 100 ab (bottom). 
Excluded regions are to the left of the solid red line. 



and 90% for RS, depending on the mass. 

To determine the reach of this analysis, we search for the local minimum of the differential 
cross section, as well as the corresponding point on the other side of the putative resonance 
peak. We then compute the integrated cross section in this region and compare it to a 
conjectured flat background. If we see a Sex excess of events, then we know to a very high 
degree of confidence that we are seeing something new, since nowhere in the SM do we expect 
to see such behavior. We present our results for the upturn analysis at 1 fb~^ at the 7 TeV 
LHC in Fig. 

The behavior is as expected for narrow resonances, where nearly all of the signal events 
are above the upturn in cross section, and so the upturn analysis works everywhere that a 
normal counting experiment would. For very broad resonances (this definition depends on 
the precise model and resonance mass, but typically involves a ratio of width to mass that is 
larger than 50% for RS and 30% for Z'), the convolution with the falling parton luminosity 
function results in fewer events above the upturn (and, in some cases, no upturn at all), and 
they cannot be reached by this analysis. 

This method is quick and easy but is overly conservative, since the background (and 
other, similar-looking types of new physics, such as contact interactions) are falling and not 
constant. Therefore, we are neglecting some broad resonances that should be distinguishable 
from background because they don't show a particularly pronounced upturn in the differential 
cross section. We expect that the LHC reach for broad resonances should be more extensive 
than that shown here. Another shortcoming of this analysis is that it fails to account for 
statistical fluctuations in determining the start and endpoints for the integration window. 
An upturn could just as well come from a downward fluctuation in the lower-invariant-mass 
region, and our analysis does not consider this possibility. 
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5.2 Maximum likelihood method 



We now present a more powerful method for discriminating resonances from contact interac- 
tions than the analysis presented above. We perform a full maximum likelihood anlysis [23] on 
the invariant mass distribution from a 1 fb~ experiment. We simulate experiments assuming 
an underlying resonance as new physics, compare the maximum likelihood ratio obtained from 
resonance and contact interaction fits to the data, and determine over what parameters we 
can distinguish the two types of physics at the 95% confidence level based on the likelihood 
ratios. We outline details of our procedure here, and we present our results in section 15. 3[ 

At very low statistics in each bin, as we expect to see if we restrict ourselves to high 
invariant mass events, we maximize the maximum likelihood function 

Hi) = Yl /(/ii, rii), (18) 

i 

where f{fi,n) = ^^/n\ is the standard Poisson distribution and the /ij are the means 
in each bin. Using the standard distribution to perform fits is not valid, as it assumes 
normally-distributed errors while the correct underlying distribution is Poisson [21]. The 
mean in each bin is found by integrating the differential cross section over the bin width. An 
unbinned analysis is also possible, but as the unbinned case can be obtained by taking the 
limit of the bin width to zero, we derive the results for a binned analysis and comment further 
on appropriate bin sizes in Appendix [Al 

The bin means are functionally dependent on the parameters of the theory; for instance, 
with contact interactions, we have /ijc = /^j(A), while for RS we have /ijr = f^i{Mg,k/Mpi), 
with the i index labeling bins. 

With few high energy events, a better statistic than is [21] 

g = -21nA = 2^(/i, -ni + riiln— ) , (19) 

where we define 

which is the maximum likelihood ratio. The statistic Q has the property that Q G [0, oo), 
analogous to the x^- Maximizing L is equivalent to minimizing Q, which is the criterion used 
for fitting, and good fits are characterized by small Q. 

There are some important differences however from x^, which are seen when we apply the 
two statistics to a single bin (take i to only have one value). Unlike the distribution, the 
mean and variance of Q is dependent on the mean of this bin [25]. For the x^, we always have 
an expectation value (x^)/d.o.f. (or (x^)/bin) of one, regardless of the bin mean. With small 
bin mean, we instead find (Q) < 1. If we now sum over many bins, the total Q statistic will 
be dependent not only on the number of bins (like the x^)? but also on the means of those 
bins. Since the Q distribution is dependent on the bin means, which in turn are functions of 
the parameters of the theory, we have a different distribution of Q for each set of parameters. 
We therefore have no choice but to generate these through Monte Carlo simulations. 
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We will assume in everything that follows that the correct underlying new physics is a 
resonance model, and we will try to distinguish it from a contact interaction. To do so, we 
compare the values of Q that come from fitting resonance (Qr) and contact interaction {Qc) 
models to any simulated dataset. We compute 

Ag = - Qr = 2 ^ (^/i/c -^i^r + nM^^. (21) 

expecting AQ < for a dataset that looks more like a contact interaction and AQ > 
for a dataset that looks more like a resonance (the subscripts c and r stand for contact and 
resonance respectively). 

Statistical fluctuations can make a contact interaction look more like a resonance and vice- 
versa, so it is not sufficient to examine the sign of AQ. Our procedure for studying a resonance 
with mass M and coupling g is as follows: 

1. Generate data D for a pseudoexperiment with 1 fb~^ at a/s = 7 TeV using a resonance 
model with parameters M and g. 

2. Perform fits of a resonance model (M^ and g^) and a contact interaction with scale 

to the data Z), minimizing and respectively. The superscript shows that this is 
with respect to the data D. 

3. Compute AQ^ = Qf - gf. 

4. To determine the values of AQ that come from a contact interaction with scale 
faking a resonance, generate data D' based on a contact interaction with scale A^. Fit 
resonance and contact parameters to D' and compute AQ'^ . Repeat many times to 
generate a probability distribution V over possible values of AQ^ • The corresopnding 
cumulative distribution, given by the integral of "P, is C. 

5. To distinguish the resonance in the data D from a contact interaction at the 95% con- 
fidence level, we require Ag^ > C^i(0.95). In words, if AQ^ is bigger than 95% of the 
AQ'^ values based on fits to an underlying conctact interaction, then we can say with 
confidence that we have a broad resonance rather than a contact interaction. 

The maximum likelihood ratio is well-known and used frequently in particle physics ex- 
periments in performing fits and determining the parameters of a theory from data, but we 
are here using it to directly compare two competing theories of new physics that may yield 
similar distributions of observables. 

5.3 Results of statistical analysis 

We now present the main results from our analysis. We apply the maximum likelihood analysis 
from the last section to the study of RS KK gravitons and non-universal B — 3L^ Z' discussed 
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Figure 9: Reliability of distinguishing B — Z' resonance at 95% confidence level from 
VV destructive (top left), LL destructive (top right), VV constructive (bottom left), and LL 
constructive contact interactions. The shaded regions, from darkest to lightest, show regions 
with reliability: > 99%, 90-99%, 70-90%, 50-70%, < 50%. Excluded regions are to the left of 
the solid red hne, while the dashed lines show cross sections: 100 fb (top), 10 fb, 5 fb, 1 fb 
(bottom). 



in section [3lj We determine over what regions of parameter space we can distinguish such 
models from contact interactions. Examining Fig. |2l it is apparent that the destructively- 
interfering contact interactions will look most like a broad resonance. Whether it is LL or VV 
that most resembles the resonance depends on the nature of the resonance model's coupling. 

The analysis outlined in section studies only a single AQbf arising from a single simula- 
tion of a resonance in 1 fb~^ of data. To quantify the effects of the statistical fluctuations over 
many different simulated experiments, we define reliability as the percentage of experiments 
in which we can distinguish a resonance from a contact interaction at the 95% confidence level. 
This essentially tells us the probability of being able to distinguish a resonance with 1 fb~^ at 
the LHC for various possible parameters. 

Our results, showing the reliability curves for distinguishing resonances from contact in- 
teractions, are presented in Fig. [10] for RS and in Fig. |9]for the non-universal Z' . Some of the 
key points from the figures that are true for both models are: 

^In this section, we do not consider the flavor-universal B — L case, as the LEP bounds imply narrow 
resonances in accessible regions of the parameter space at the LHC. 
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Figure 10: Reliability of distinguishing RS resonance at 95% confidence level from VV destruc- 
tive (top left), LL destructive (top right), VV constructive (bottom left), and LL constructive 
contact interactions. The shaded regions, from darkest to lightest, show regions with reliabil- 
ity: > 99%, 90-99%, 70-90%, 50-70%, < 50%. Excluded regions are to the left of the sohd red 
line, while the dashed lines show cross sections: 100 fb (top), 10 fb, 5 fb, 1 fb (bottom). 



1. It is most difficult to distinguish a resonance from a contact interaction that interferes 
destructively with the SM (in agreement with our earlier reasoning). For Z', it is the 
VV destructive contact interaction that is most likely to fake the resonance since that 
is the nature of the Z"s own interference and coupling. For RS, where SM interference 
does not modify the total cross section, both destructive models do equally well, with 
only modest differences between the others. 

2. For narrow resonances {k/M-p\ < 0.4, e < 0.75), the reliability curves track the cross 
section contours. Here, the relevant question is only the number of signal events required 
for detection: with five or more events, we can reliably distinguish between a resonance 
and contact interaction. 

3. For broad, low-mass resonances {k/Mpi > 0.4 and M < 2 TeV for RS, e > 0.75 and 
M < 2.5 TeV for Z'), we can universally discriminate between contact interactions and 
resonances because of the large number of events at 1 fb~^. The effects of statistical 
fluctuations become negligible with many events. 
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Figure 11: Plots showing mean best-fit values of A for RS models compared to VV (left) and 
LL (right) destructive contact interactions. Contours of A are shown in solid orange lines. 
They are for the left plot: 7 (top), 8.5, 10, 11.5 TeV (bottom); for the right plot: 5 (top), 6, 
7, 8 TeV (bottom). The reliability regions are shaded and cross sections shown with dashed 
lines, as in Fig. [101 



4. As expected, broad, high-mass resonances (/c/Mpi > 0.4 and M > 2 TeV for RS, e > 0.75 
and M > 2.5 TeV for Z') are the hardest to distinguish from contact interactions at lower 
statistics. We find that it is still possible to detect broad resonances with confidence 
over most of the kinematically accessible parameter space, but the reliability of doing so 
diminishes when the width is greater than F/M ~ 20%. 

The analysis described above was done only for a single, isolated resonance. Generally, 
when we encounter strongly-coupled physics, we may find multiple resonances. We certainly 
expect to see additional resonances in the RS framework, namely the higher KK modes. The 
interference between resonances may be sufficiently high that we can no longer neglect the 
heavier modes, depending on resonance width and spacing. 

As a check, we compute the contribution of multiple graviton KK modes (and their mutual 
interference) to the Drell-Yan process in RS. For k/Mp\ < 0.5, the effects of higher KK modes 
on the differential and total cross sections around Mg are minimal. At higher coupling, we see 
enhancements to the cross section of up to ~ 15 %; most importantly, however, the shape is not 
distorted. Thus, these effects do not substantially change the outcome of the single-resonance 
analysis; the marginal increase in event rate will make at best make model discrimination 
slightly easier. 

To better understand how the fitting procedure works, we plot the mean best-fit values of 
A for RS resonances in Fig. [TTl both for LL and VV destructive contact interactions. It is 
apparent that the best-fit values of A are strongly correlated with the total number of events. 
The contact interaction most likely to fake a broad resonance is one that predicts the same 
number of events as the resonance, only with the events distributed in a different way. This is in 
accordance with our picture of unfortuitously distributed upward- and downward-fluctuations 
making a contact interaction look like a resonance. 

In our analysis, we choose a bin size of 100 GeV, since that is on the order of the detector 
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Figure 12: Plots showing dimuon cross section for the first KK mode of the graviton in RS 
models (to muons, electrons, and photons) at y/s = 7 TeV. Legend: green is 1 pb, blue is 100 
fb, orange is 10 fb, purple is 1 fb. The black curve indicates the cross section for 5 events. 



resolution for dimuon invariant masses around 1 TeV [26]. The range of invariant masses 
that we consider is 300 GeV to 4 TeV|f| We have not included detector smearing effects. A 
discussion of the effects of including detector smearing and varying the bin size (justifying a 
binned analysis) is found in the Appendices. 

5.4 Electrons and photons 

The non-universal Z' we have been considering evades the LEP bounds by coupling only to 
muons, and there will be no excess dielectron or diphoton final state events over the SM 
prediction. RS gravitons, however, couple to both electrons and photons, and these should 
make it easier to distinguish gravitons from contact interactions. 

Electron and photon final states, like muons, have small SM backgrounds and are straight- 
forward to detect. Depending on cuts imposed to minimize fake rates, detector efficiencies for 
electrons and photons can be as high as 60-85% [H] . The cross section for photon production 
by a KK graviton is twice that for muons, and so the total event rate increases by just under 
a factor of four when we include all channels (it is smaller than four because of the lower 
efficiencies for photons and electrons). 

We plot our results for the RS cross section to electrons, muons, and photons in Fig. [121 
We do not repeat the full statistical analysis and comparison to contact interactions because a) 
we expect the qualitative result to be the same, and b) it is difficult to systematically explore 

^Since the area around the Z peak is used to normahze the cross section and match the PDFs, we choose 
a lower value that safely avoids this region 
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the full range of possible contact interactions, since we could have different coefficients for the 
operators involving different final states, etc. However, we do compare the cross section to 
that for muons alone and find that the new cross section contours are now shifted upward 
in mass by about 400 GeV. We therefore expect to be able to probe resonances to 400 GeV 
higher masses than if we relied only on muons, and the reliability regions would also shift 
accordingly, meaning that we have difficulty discriminating 2.8 TeV gravitons from contact 
interactions, rather than the 2.4 TeV gravitons we had with just muons. 

6 Angular distribution 

In this section we demonstrate how to use the pseudorapidity (77) distribution of the final state 
high Prp muons to provide another handle in discriminating between models. We define a new 
variable called ellipticity, which measures how longitudinal the outgoing muons are in either the 
forward or backward directions, with no need to explicitly identify the direction of the quark 
vs. antiquark (as would be necessary for forward-backward asymmetry determination). We 
will see that this is a more efficient way to identify parity violation in new physics, particularly 
for broad resonances. 

The SM violates parity, whereas new physics models may not. In a pp collider, this is 
encoded in the shape of the pseudorapidity distribution, as can be seen in Fig. [131 which 
shows the hadronic differential cross sections {da/dr]'^) for the standard model, where we have 
defined 77^ as the pseudorapidity of the outgoing fi"^ from pp — The outgoing /i~s are 
preferentially scattered forward relative to the quark in the SM, and since the quark could 
come from either proton, the r)~ distribution for the SM demonstrates a characteristic double 
peak structure. On resonance, new physics contributions that respects parity, such as the RS 
model, would wash out this effect, resulting in a transversely peaked distribution as shown on 
the right side of Fig.[T31 Note that, at a pp collider, studying the forward backward asymmetry 
in the CM frame will not capture this effect since, without identifying the quark direction, the 
symmetry of the 77^ distribution would make the distinguishing peak structure cancel. 

The shape of the 77" distribution is significantly different than that of the t/"*" distribution 
and requires explanation. The underlying hard interaction involves a valence quark whose 
momentum fraction is generally much larger than the antiquark's and the resulting collision 
is boosted in the quark's direction. Furthermore, the Standard Model preferentially scatters 
fi~s into the same direction as the quarks as shown in Fig. [TU and when combined with 
the effect of the large boost, the resulting 7/~ distribution is peaked away from the zero as 
shown by the dashed lines in Fig. [131 Since the valence quark could come from either proton, 
the distribution must be symmetric in 7/~ and the two dashed curves are added to give the 
final double arched curve in Fig. [121 The /i"*" cross section is more centrally peaked since the 
H'^ is preferentially scattered backwards from the quark and, after the boost, becomes more 
transverse with smaller rapidity. A pp collider would separately show the effects of parity 
violation in both of the rj distributions since these curves would are symmetric; however, the 
large boost at a pp collider magnifies the effect in the 77" distribution while simultaneously 
diminishing the effect in t/"*". 
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Therefore, we study a new observable, the ellipticity (Erj), to probe the shape of the ri~ 
distributions and we demonstrate its usefulness in discriminating between new physics models 
at the LHC. We calculate the ellipticity with the formula 
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where r]^ is the pseudorapidity of the outgoing /x^. This quantity seeks to exploit the distinc- 
tive shape of the ri~ distributions. The end result is relatively insensitive to values of a; ~ 1 
and ?7inax ~ 2.5, and so we have chosen x = 1.0, r/max = 2.5 for this study. Ultimately, the 
value of X could be optimized. Furthermore, studying a ratio has the usual advantage that 
some systematic uncertainties, such as the K-factor, integrated luminosity, PDF uncertainties, 
etc., will be common to both quantities and thus divide out. 

We will also contrast ellipticity with the center-edge asymmetry [29] (Ace), a related 
quantity based on the angular distribution in the CM frame that has been considered for 
spin determination. The center-edge asymmetry is defined using the quantity z = cos 6* = 
tanh( ^ "'^ ) and is 
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This quantity attempts to exploit the differences in the cos^* distribution due the spin of 
the structure of the signal. For our study, we have taken z^^^ = 1.0 and 2;* = 1/2. A more 
thorough study of the variation of z* was carried out in [29]. Also, a modified version of 
center-edge asymmetry (Acc) was proposed in [30] that is based on A?7 = ri~ — j]^ instead of 
cos^*. Since cos 6** = tanh(A?7/2), A^e and A^^ behave similarly and any conclusions we draw 
for will also apply to Ace. Notice that while Ace and are defined analogously, with 
z ri~ . We will demonstrate that, while E^ and Ace are both sensitive to the underlying 
spin structure, Er, provides more discriminating power among models with different parity 
violating character. 

Spin determination studies of narrow resonances have focused on the distribution of cos 9* = 
tanh(?7~ — t]~^)/2, where 6* is the center of mass frame partonic scattering angle between the 
quark and the n~ . Quantities that rely on cos^* differ from quantities that rely on 77" for 
several reasons. 

• cos^* depends on a difference of rapidities, which makes the quantity boost invariant. 
The boost magnifies the effects of parity violation and, since cos 6* is boost invariant, the 
characteristic double peak feature of the standard model 77" distribution is not present 
for cos 6* . 

• The forward-backward asymmetry {Afb) depends on cos 6'* and requires that one deter- 
mine the direction from which the quark originated. This direction is not known in a pp 
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Figure 13: The left plot shows the hadronic differential cross section in r/" (upper) and rj^ 
(lower) for pp — )■ fj^~^fi~ via the standard model at Ecm = 7 TeV with M^^j > 400 GeV. The 
dashed curves represent the individual contribution from each of the two terms in Eq. [21] 
and the solid curve is the sum. The right plots shows the same cross sections for the SM 
(sohd black), RS model with Mg = 1300 GeV, k/Mpi = 0.3 (solid, red), {B - 3L^) Z' model 
with Mg = 1300 GeV, k/Mpi = 0.48 (dashed, green), LL composite model with destructive 
interference and A = 4550 GeV (dotted,blue), and VV composite model with destructive 
interference and A = 6150 GeV (dot-dashed, orange). They have been normalized to compare 
shape. 



collider a priori and the assignment of 6* for each event is ambiguous. For very forward 
events, the boost required to bring the muons to the CM frame aligns with the direction 
of the quark's momentum [28]. A rapidity cut to remove the low rapidity events would 
generate a sample for which we could measure ApB] however, it would also remove a 
majority of the signal events, which tend to be centrally peaked, and thus would not be 
useful until a significant amount of data has been recorded. 

Another version of ^4^?^ was proposed in [30] which probes differences in the distribution 
of |?7~| — |?7'''|. Since this variable is sensitive to the boost of the system, it behaves 
similarly to E^. Moreover, they showed that the inclusion of transverse events did not 
significantly affect the discriminating power of this variable and that it does not suffer 
the deficiency of the traditionally defined version of ApB based on cos^*. 

• Acc cannot distinguish between models with the same spin structure, but with different 
amounts of parity violation. For example, the composite models being considered arise 



22 



Figure 14: The normalized differential partonic differential cross section for uu — > jj,~^fi~ 
via the standard model is shown in both plots (solid black). The left plot shows the cross 
section for the SM with RS graviton of mass Mg = 1300 GeV and coupling fc/Mpi = 0.3 for 
the processes uu — )■ (dotted red), and gg G* ^ /^^/^~ (dot-dashed purple). The 

right plot shows the differential cross section for the process uu — )■ for the SM plus 

{B - 3L^) Z' with mass Mz' = 1300 GeV and coupling e = 0.48 (dotted green), SM plus 
LL (destructive) composite model with A = 4550 GeV (dashed orange), and SM plus VV 
(destructive) composite model with A = 6150 GeV (dashed orange). All plots are shown with 
partonic CM energy a/I = 1100 GeV. 



from integrating out a heavy spin one resonance, and so they share the same initial state 
and spin structure to the Z' models. After convolving the partonic differential cross 
sections with the parton luminosity function, any parity violating effects are washed out 
as demonstrated in Fig. [151 These models are indistinguishable using Ace- Conversely, 
not only can E^^ can distinguish models with similar spin structure, but with different 
amounts of parity violation. 

• A difference in Acc indicates either a different spin structure and/or coupling to initial 
state gluons. Thus Ace can identify resonances mediated by an RS graviton, but it 
cannot distinguish between the rest of the models that we considered. In this regard. 
Ace and E^ are similar. 

The calculation of the ellipticity is straightforward. The following expression gives the 
differential cross section in M^^ and y^, where y^ is the rapidity of fi^: 

da 1 1 1 J _y ^ dcxij 





(24) 



We have used the definitions 
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Figure 15: . The (hadronic) differential cross section in cos 6'* for pp — )■ fi^fi^ is shown at 
Ecm = 7 TeV with M^^ > 400 (left) for the same models as Fig. [HI The right plots shows 
the same cross section for the same models as Fig. [T7| at Ecm = 10 TeV with M^^ > 1000. 
The RS model corresponds to the dashed red curve and shows a distinctive shape due to the 
spin-2 character of the graviton and its coupling to initial state gluons. The other models are 
shown in solid black. 



and the allowed region is 

-log^<,- + ,-<log^ 

— oo < y~ — < CO. (26) 

The two terms in braces in Eq. [2^ represent the quark coming from either proton. The muon's 
rapidity can be replaced with its pseudorapidity {y"^ — ^ ^7^) at the energies we are considering, 
Inserting Eq. flM|) into Eq. fl22l) gives the average value of the ellipticity. We similarly calculate 
the center-edge asymmetry using 

d'^a 1 -^--^ dCijir) daijirs) 



r = M^IEL, (27) 



dW^dcosO* E}^'^ dr dcosO 

which is then inserted to Eq. f l23|) . Finally, when integrating over the invariant mass M, 
one chooses different integration regions depending on the width of the resonance. When the 
resonance is narrow, events are selected in a single invariant mass bin around the resonances, 
which enhances signal over background. However, for broad resonances, there isn't a conve- 
nient mass window that enhances the sample, and one is forced to include all of the events 
above a given invariant mass, M^^, which will vary according to the shape and mass of the 
resonance. We have chosen the values of M,™* to be 400 and 1000 GeV for the 7 and 10 TeV 
runs, respectively. 

The appearance of (parity conserving) new physics above the mass cut will give larger 
values for the ellipticity due to the larger proportion of of centrally peaked fi~s. For example, 
the representative RS model has a contribution to the DY cross section from the initial state 
gluons. These gluons tend to have comparable momentum fractions and so there isn't a large 
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Figure 16: A^c (left) and E^^ are shown for 1000 pseudo-experiments for Ecm = 7 TeV, J dtC = 
1 fb-i with M^^ > 400 GeV. The models shown are a RS model with Mg = 1300 GeV, 
k/Mpi = 0.3 (red), {B-3L^) Z' model with Mg = 1300 GeV, e = 0.48 (green), LL composite 
model with destructive interference and A = 4550 GeV (blue), and VV composite model with 
destructive interference and A = 6150 GeV (orange). 



boost as in the qq initial state; therefore, the ri~ distribution from gg tends to be centrally 
peaked. The contribution from qq ^ G* ^ A*^/^" shows a slight asymmetry, but it is dwarfed 
by the gluon contribution that washes out any distinctive shape. 

The rather high invariant mass cut means there are only 100 or so events for 1 fb~^ at 7 
TeV and so the observables will have large statistical uncertainties. To model this uncertainty, 
a distribution of Ace and Erj was generated for 1000 pseudo-experiments for 1, 5 and 10 fb~^ at 
various CM energies. A representative set of distributions are shown in Fig. [161 The pseudo- 
experiments were generated with Pythia 8 [331 El]- The LHC will only be able to make a 
single measurement and the distributions of Fig. [16] for the completing models are too similar 
for effective discrimination. There is some separation of the RS model from the rest, but to 
say the models are distinguishable is a bit of a stretch since there is significant overlap with 
the other distributions. We conclude that for masses greater than 1000 GeV, 1 fb~^ at 7 TeV 
is not enough integrated luminosity to be able to reliably distinguish among the models using 
these observables. Indeed, it is unlikely that any angular distribution observable will be able 
to discriminate models with so few statistics. 

Discriminating power improves with more events, whether due to higher cross section (due 
to lighter new physics or higher CM energy) or higher integrated luminosity. The distributions 
are presented in Fig. [T7] for E^m = 10 TeV with 1 fb~^ (upper plots) and 5 fb~^ (lower plots) of 
luminosity for the same models as in Fig. [161 but with modified couplings and higher resonance 
masses or composite scales. The distributions in the lower plots show considerable separation 
indicating that for certain parameters, these models are distinguishable. For example, if the 
experiment measured a value of E^ of 0.5, there is little chance that the LL composite model 
could be the underlying physics. The Acc distributions shows clear separation between the 
RS model and the other spin 1 models, whereas the E^^ distribution distinguishes between the 
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Figure 17: Ace (left) and -E"^ are shown for 1000 pseudo-experiments for Ecm = 10 TeV, 
/ dtC = IJb-^ (upper) and 5 fb'^ (lower) with M^^ > 1000 GeV. RS model with Mg = 2000 

= 0.61 (green), LL 
and VV composite 



GeV, k/Mpy = 0.5 (red), {B - 3L^) Z' model with Mg = 2000 GeV, e 
composite model with destructive interference and A = 5940 GeV (blue), 
model with destructive interference and A = 8100 GeV (orange). 



RS model, the LL composite model, and the parity conserving Z' and VV composite models. 
This is a clear demonstration that E^j has more sensitivity to parity violation than A^c and will 
be a useful observable for model discrimination for broad resonances. Finally, Acc and E^^ are 
complementary variables. Measuring A^e will help discriminate between models of differing 
spin structure, after which a measurement of E^ will further differentiate models based on their 
parity violating character. Once the models are distinguished, a more thorough investigation 
of the properties of the resonance can be performed. 



7 Conclusions 

In gearing up to study new physics at the LHC, we must consider the different features of the 
new collider. We have shown that, even in the first LHC run, its energy advantage over the 
Tevatron will allow it to probe new resonance physics, but only when couplings are sufficiently 
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large to give sufficiently many signal events at relatively low luminosity. This has led us to 
the search for broad resonances. We have seen that careful statistical analyses will let us not 
only search for new resonances, but also discriminate them from contact interactions (as well 
as background). 

We have furthermore considered angular distributions, a discriminator most likely to be 
useful only when the LHC runs at higher energy and luminosity. We demonstrated that the 
ellipticity is sensitive to parity-violation, even without knowing the quark direction, whereas 
the center-edge asymmetry is only sensitive to the spin structure of the interaction. Both 
observables provide complementary, useful, information. 
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Appendices 

A Binning 

We performed a binned analysis in section 15. 3[ which introduces uncertainties relative to an 
unbinned analysis that are proportional to the square of the bin size [27]. Nevertheless, we 
anticipate that the binned analysis gives comparable results to an unbinned analysis for the 
following reasons: 
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Figure 18: Comparison of statistical analysis with different bin sizes for a Z' compared to a 
VV destructive contact interaction. The reliability shadings and cross sections are as in Fig. [9], 
and the analysis was performed with bins of 50 GeV (left) and 25 GeV (right) width. 
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1. The features of a broad resonance are on a scale ~ F, which is must larger than the 100 
GeV bin width. 

2. With few events at high invariant mass, we never expect to see more than a single event 
(or two at most) per bin. 

3. At low invariant mass, where we expect more events, the distribution is dominated by 
SM interference, which we expect to be similar for different competing models of new 
physics (such as as Z' vs. a destructively-interfering contact interaction. A counting 
experiment suffices in this region, as we don't need high sensitivity to the shape of the 
invariant mass distribution here. 

We present the results of a binned analysis with bin sizes 50 GeV and 25 GeV in Fig. [181 
which are reduced relative to our standard 100 GeV bin size. It is evident that any differences 
between them are minor and do not affect our results qualitatively; therefore, the 100 GeV 
bin size is sufficiently small to give accurate results. 

Reducing the bin size does have some advantages, however. In particular, the accuracy of 
the fit parameters is enhanced by performing an unbinned analysis, even though AQ does not 
change much between them. In our case, however, we are not actually concerned with a precise 
determination of the resonance parameters, but whether characteristics of the invariant mass 
distribution allow us to generically distinguish it from other types of new physics, and this 
can be achieved accurately using a binned analysis. 

B Detector smearing 

The muon detectors at ATLAS and CMS have some finite resolution that introduces uncer- 
tainties into the invariant mass measurements. This has the tendency to smear the resulting 
distribution that comes out of the detector. The smearing should be applied both to the 
differential cross section used to generate the simulated data, as well as the contact and res- 
onance distributions used to fit the data. This makes the whole simulation/fitting procedure 
much more computationally-intensive, particularly considering the number of fits that are re- 
quired to scan over the parameter space and generate the AQ distributions from the contact 
interactions. 

Smearing tends to give otherwise narrow resonances a width of about the detector reso- 
lution, which at (9(100 GeV) is still too small to leave any doubt that the new physics is a 
resonance. For broad resonances, however, it could further obscure the existing features of the 
resonance and make it look even more like a contact interaction. To determine the effects on 
our result, we repeated the analysis for a Z' assuming Gaussian smearing using the detector 
resolution given by [26] . Our result is shown in Fig. [T9l 

The introduction of detector smearing does change the results, albeit by a relatively small 
margin (for broad resonances, we have 50 GeV loss in discriminatory power at fixed e). The 
detector effects do not, however, qualitatively change the results in any way, introducing only 
a small, horizontal shift in the reliability curves. Given that a proper analysis at ATLAS or 
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Figure 19: Comparison of statistical analysis without (left) and with (right) Gaussian smear- 
ing of the final muon states for the non-universal Z' compared to a VV destructive contact 
interaction. The reliability shadings are as in Fig. |9l as are the cross sections. 



CMS will require full-scale modeling of the detector, we are content to accept the non-smeared 
results as valid to within such a ~ 50 GeV window. 
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